Magneto-optics of bilayer graphene inclusions in rotational-stacked 

multilayer epitaxial graphene 
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Additional component in multi-layer epitaxial graphene grown on the C-terminated surface of SiC, 
which exhibits the characteristic electronic properties of a AB-stacked graphene bilayer, is identified 
in magneto-optical response of this material. We show that these inclusions represent a well-defined 
platform for accurate magneto-spectroscopy of unperturbed graphene bilayers. 

PACS numbers: 71.55.Gs, 76.40.+b, 71.70.Di, 78.20.Ls 



I. INTRODUCTION 

The absence of the energetic gap in the excitation spec- 
trum of graphene^ is considered as a possible drawback 
preventing the straightforward application of this emerg- 
ing material in electronics. This is despite numerous 
efforts, such as those implying surface patterning^ and 
substrate- or adsorbents-induced interactions.— ~— The 
possibility to open and tune the band gap in the bilayer 
graphene has recently been demonstrated by applying an 
electric field perpendicular to the graphitic planes 6 - - — and 
this is a key element to construct a transistor, the build- 
ing block of electronic circuits. The band gap engineering 
is "typical" of the bilayer and is not reported in tri- and 
more-layer graphene specimens where semi-metallic be- 
havior dominates^ From the viewpoint of applications, 
the bilayer graphene thus becomes almost equally appeal- 
ing material as graphene itself. 

Optical spectroscopy has played an important role in 



resonance signal of a small flake using the gate-controlled 
differential technique. The optical response at a fixed 
magnetic field was then studied as a function of the car- 
rier density. Such differential spectroscopy was efficient 
in case of exfoliated graphene monolayers 1 1 i 18 but it pro- 
vides more complex results when applied to the bilayer 
graphene. In this latter system, the change of the gate 
voltage affects not only the carrier density but also mod- 
ifies significantly the band structure and data interpreta- 
tion is by far more elaborated. 20 

In this paper, we demonstrate that certain class of 



previously reported AB-stacking faults 
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stance, it allows to directly visualize the electric-field in- 
duced energy gap in this systemi 7 -^ On the other hand, 
only relatively scarce information has been up to now col- 
lected from magneto-optical measurements^ 2 - This fact 
might be surprising when noticing the potential of Lan- 
dau level (LL) spectroscopy which has been widely ap- 
plied to other graphene-like systems. Magneto-optical 
methods have, for example, convincingly illustrated the 
unconventional LL spectrum in graphene, have offered a 
reliable estimate of the Fermi velocity or invoked the spe- 
cific effects of many-body interactions between massless 
Dirac fermionsJ^— 

So far, the only magneto-optical experiments on the 
bilayer graphene have been reported by Henriksen et 
al.r^ who succeeded to probe a relatively weak cyclotron- 



wise rotationally-ordered multilayer epitaxial graphene 
(MEG) / 3 ' 25 — show the characteristic features of well- 
defined graphene bilayers. These inclusions, identi- 
fied here in magneto-transmission experiments, repre- 
sent therefore a suitable system for accurate magneto- 
spectroscopy studies of unperturbed bilayer graphene. 



II. SAMPLE PREPARATION AND 
EXPERIMENTAL DETAILS 

The growth of MEG samples studied here was per- 
formed with a commercially available horizontal chem- 
ical vapor deposition hot-wall reactor (Epigress V508), 
inductively heated by a RF generator. Epitaxial MEG 
films were grown on semi-insulating 4H-SiC(0001) on- 
axis C-terminated substrates at 1600°C in Ar atmo- 
sphere. The growth rate was controlled by the Ar pres- 
sure (~100 mbar) which was found to directly influence 
the evaporation rate of Si atoms. 

To measure the infrared transmittance of our samples, 
we used the radiation of a globar, which was analyzed 
by a Fourier transform spectrometer and delivered to the 
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FIG. 1: (color online) Transmission spec- 
tra of MEG with ~100 layers recorded 
at the selected magnetic fields below 12 
and above 16 T in parts (a) and (b), re- 
spectively, in both cases above the rest- 
strahlenband of SiC Whereas the transi- 
tions marked B to I correspond to electri- 
cally isolated graphene monolayers, tran- 
sitions denoted n = 1 to 5 match to 
inter-LL transitions in an unperturbed 
graphene bilayer, following the coding 
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cessive spectra in parts (a) and (b) are 
shifted vertically by 0.14 and 0.23, re- 
spectively. 



sample via light-pipe optics. The transmitted light was 
detected by a composite bolometer kept at T = 2 K and 
placed directly below the sample. Measurements were 
carried out in superconducting (B = — 13 T) and resis- 
tive {B = 13 — 32 T) solenoids with spectral resolution 
of 0.5 and 1 meV in the range of magnetic field below 
and above B = 13 T, respectively. All presented spectra 
were normalized by the sample transmission at B = 0. 

The samples were characterized in micro- Raman scat- 
tering experiments which, similarly to previous studies,— 
revealed, depending on location, single-component 2D 
band features, characteristic of graphene simple elec- 
tronic bands and of decoupled graphitic planes in multi- 
layer epitaxial graphene grown on the C-face of a SiC sub- 
strate, or multi-component 2D band features characteris- 
tic of Bernal stacked graphite. In this paper, we present 
transmission data obtained on one particular specimen 
with a high number of graphitic layers (~ 100) and grown 
on a SiC substrate with a reduced thickness of ~ 100 /im. 
Due to this latter condition, the spectral region of total 
opacity of the sample only covers the SiC reststrahlen- 
band (~85-120 meV), i.e., it is significantly narrower 
as compared to the case of the preceding studies . 28 ' 29 
In spite of these efforts to expand the available spec- 
tral range, a relatively weak transmission was still found 
around the energy of 200 meV, due to double-phonon ab- 
sorptions in the underlying SiC substrate and transmis- 
sion spectra are affected by strong interference patterns 
due to the relatively thin substrate. These two effects 
prevented measurements in the energy range below the 
reststrahlenband of SiC. 



We denote those lines by Roman letters, following the 
initial work and notation of Sadowski et al*££& These 
dominant spectral features are equivalent to the charac- 
teristic lines observed in the magneto-optical response of 
exfoliated graphene monolayers! 14 i 15 i 18 The subsequent 
absorption lines labelled here as B — > I correspond to 
transitions L_ m /_ m _i) — s-L m+1( - m ) with m = 0->7 be- 
tween LLs in graphene: E rn — s\gn(m)E\y/\m\, where 
Ei = vpy/2h\eB\. The apparent Fermi velocity is ex- 
tracted to be v F = (1.02 ± 0.02) x 10 6 m.s -1 . Intrigu- 
ingly, the L-i(o) — >Lrjm transition exhibits a significant 
broadening above 16 T, which could be tentatively re- 
lated to electron-phonon interaction. This effect will be 
discussed elsewhere. 

The main focus of the present work are other spec- 
tral features, i.e., the transmission dips denoted by 
n=l, n= 2, . . . n—5 in Fig. Q] These absorption lines 
are significantly weaker than the dominant "graphene 
lines", but are still well resolved in our spectra. As it 
can be seen in Fig. [5J in contrast to the dominant tran- 
sitions, these weaker lines follow a nearly linear in B 
dependence. As this behavior is characteristic of massive 
particles and because graphene bilayer is the simplest 
graphene based system with such particles, we antici- 
pate that electronic excitations within graphene bilayer 
inclusions are responsible for the n=l, n=2, . . . n=5 
transitions. The energy ladder e n ^ of LLs in a graphene 
bilayer can be easily calculate d 30 ' 31 within the standard 
four band model which only considers the two most rel- 
evant coupling constants 70 and 71 (se e.g. Ref. |32| for 
their definitions): 



III. RESULTS AND DISCUSSION 

Typical transmission spectra of the investigated sam- 
ple are shown in Fig. [TJ The dominant absorption lines 
which are observed in these spectra show the characteris- 
tic \/i3-dependence (see Fig. [5]) and correspond to inter- 
LL transitions in electrically isolated graphene sheets. 
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Here, a positive (negative) integer n indexes the elec- 
tron (hole) LLs. \x = —1 accounts for the topmost 
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valence- and the lowest conduction-band, whereas fi = 1 
corresponds to two other, split-off bands. As illustrated 
in the inset of Fig. [21 optically active inter Landau level 
transitions in a graphene bilayer fulfill the \n\ — > \n\ ± 1 
selection rule. The energies of such transitions are plot- 
ted in Fig. [2] with black solid lines. Those lines account 
for the transitions within the fj, = — 1 bands. To re- 
produce the experimental data, we have adjusted the 71 
parameter whereas the Fermi velocity vp which defines 
Ei (i.e. the intra-layer coupling 70 = 3150 meV) has 
been fixed at the value derived from the monolayer-like 
transitions. A fair agreement is obtained between the 
calculated (solid lines in Fig. [2j) and measured energies 
of n=l, n=2, . . . n=5 transitions. Optical absorptions 
involving LLs of higher indexes (e.g. n=6, n=7 and 
n=8, see Fig. [5]) could also be observed in the spectra, 
nevertheless, these lines are very weak and only visible 
in a limited range of magnetic fields. Traces of inter-LL 
transitions due to split-off bands (/1 = 1 in Eq. [T]), can 
be also identified in our data and experiments focused on 
this particular set of transitions are in progress. 

A pronounced departure of the observed bilayer tran- 
sitions from the linearity in B clearly shows the limits 
of the parabolic approximation which is often used for 
graphene bilayers in the close vicinity of the k — point, 
and within which the LLs are strictly linear with the 
magnetic fiek h 32 ' 33 As can be seen in Fig. [21 the posi- 
tions of all these lines can be very well reproduced with 
a parameter 71 — (385 ± 5) meV, and these experiments 
thus refine the value of this parameter reported previ- 
ously from optical studies at zero magnetic field^ ' 10 ' 11 
The intriguing splitting of the n=l and of the n=2 lines 
at high magnetic fields is beyond our simple model and 
will be discussed later on. 

The simplified model of LLs in the pristine bilayer 
graphene provides us with reasonably accurate descrip- 
tion of its magneto-optical response, even though it ne- 
glects the electron-hole asymmetry (mainly induced by 
tight-binding 74, A' parameters) j 6 ' 8 ' 10 ' 11 as well as a pos- 
sible gap opening at the charge neutrality point. Never- 
theless, it should be noted that the optical response of 
the graphene bilayer has only been unambiguously iden- 
tified above the reststrahlenband of the SiC substrate 
and therefore, we cannot exclude a possible appearance 
of an energy gap, up to a few tens of meV, at the k = 
point. For the same reason, we can only estimate a very 
higher limit for the carrier density in the studied bilayer 
of 2 x 10 12 cm -2 . However, the real carrier density is very 
likely similar to that of the surrounding (electrically iso- 
lated) graphene sheets, i.e. below 10 10 cm -2 , as reported 
in Refs. [TH andl22l. We also point out that relatively nar- 
row linewidths of the order of 10 meV (relaxation time 
in sub-picosecond range) serve as an indication of rather 
high electronic quality of these bilayer inclusions, com- 
parable or even better than other bilayer systems . 12 ' 34 ! 35 

Equivalent bilayer-like spectral features are recurrently 
identified practically in all studied specimens, neverthe- 
less, with a strongly varying intensity. In general, we 



can say that the relative intensity of these "bilayer" 
lines increases with the total number of layers in MEG 
and these transitions are practically invisible in spec- 
imens with less than 10 layers reported in very first 
magneto-spectroscopy studies! 13 ' 28 This finding serves 
as an indication that we indeed observe graphene bi- 
layer inclusions and not regions of a local AB-stacking 
which might be also speculated to appear in rotation- 
ally stacked multilayers. Such Moire-patterned AB- 
stacked areas have been recently visualized in MEG by 
STM/STS measurements^.^ 8 - We further assume that 
twisted graphene layers which results in the Moire pat- 
terned bilayer should not provide us with so well-defined 
AB stacked bilayers as we observe in our data. Let us 
also note that if we compare the relative intensity of ob- 
served transitions, we can roughly estimate that in none 
of the investigated samples the ratio between bilayers and 
monolayers exceeded 10%. 

We should also emphasize that the appearance of 
Bernal-stacked faults in MEG, which have a form of well- 
defined bilayers, is not a signature of bulk graphite. In 
this well known material, the if-point electrons indeed 
mimic massive carriers in the graphene bilayer, but with 
an effective inter-layer coupling 271 instead of 71 in a 
real graphene bilayer] 31 ' 39- — This twofold coupling in the 
effective bilayer model for K point electrons implies a 
characteristic effective mass twice enhanced in compari- 
son to that of massive Dirac fermions in true graphene 
bilayer and consequently, also a twice lower energy sep- 
aration between adjacent interband inter-LL transitions, 
cf. Fig. [2] of this paper with the fan chart in Ref. |40T 

The remaining unclarified point of our study is the 
splitting of the bilayer lines, which is clearly visible for 
transitions n=l and n=2 around B — 17 and 26 T, 
respectively. In the following, we discuss two different 
scenarios for this splitting. One possible explanation in- 
vokes the electron-hole asymmetry, reported recently in 
graphene bilayers graphene! 6 ' 10 ' 11 Based on this assump- 
tion, the magnitude of the splitting for the n-th transi- 
tion, relative to the transition energy is expressed bys^ 2 - 

2(A77i + 274/70) 
yjn{n + l) + v /(n + l)(n + 2)' 

For the values A' = 0.02 eV, 71 = 0.4 eV, 74/70 = 
O.OSj 6 -"— i 10 ' 11 our measured value for n = 1 (about 0.08) 
is very well reproduced. However, the splitting due to 
electron-hole asymmetry should be seen for all magnetic 
fields, while, as can be seen in Fig. [TJd, we only observe 
it in a relatively narrow range of B. 

Perhaps a more natural explanation for this line split- 
ting would be an avoided crossing between the transition 
L-n(-n-i) L n+ i( n ) and some other transition with a 
much smaller oscillator strength, so that it is not seen 
far from the crossing point. One can see directly from 
Fig. [2] that the bright transitions n = 1,2 are crossed 
by the dark (i.e., dipole-inactive in case of a zero trigo- 
nal warping) transitions L (_4) — > L 4 ( ), L 0( -_ 7 ) — > L 7 ( ), 
respectively, approximately at the observed values of B 
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FIG. 2: (color online) Fan chart: Points marked with Roman 
letters, having a characteristic \/B-dependence, correspond to 
inter-LL transitions in electrically isolated graphene sheets^ 
Points denoted by index n represent inter-band inter-LL 
transitions in the graphene bilayer L_ n (_„_i) — >L n+1 ( n ), as 
schematically shown in the inset. The full gray lines show 
expected energies of transitions in the graphene monolayer 
for vp — 1.02 x 10 6 m.s -1 , full black lines correspond to pre- 
dicted positions of the absorption lines in the graphene bilayer 
(only parameters 70 = 3150 meV and 71 = 385 meV consid- 
ered). The dashed lines denote theoretical positions of two 
trigonal-warping-induced transitions in the graphene bilayer 
L (_4) — >L 4 ( ) and L (_7) — >L 7 ( ). 

(the crossing occurs at a very sharp angle, which brings 
a significant uncertainty). These transitions are allowed 
only due to the presence of the trigonal warping of the 
electronic bands, which mixes levels L m with \m\ differ- 
ing by an integer multiple of 3, see Ref. [3(i The ratio of 
the oscillator strength of the L (_4) — > L 4 ( ) transition 
to that of the bright n = 1 transition can be estimated 42 

as (25/108)(737i/7o) 2 (VftvF) 2 - °- 02 at B = 25 T - 
The L (_4) —¥ L 4 ( ) transition is therefore not expected 
to be seen in the experiment unless some other, pos- 
sibly resonant, admixture mechanism is taken into ac- 
count. Coupling between L (_4) — > L^q) and n = 1 
transitions should be quite strong as the observed "anti- 
crossing splitting" is of about 20 meV. 



We have speculated this mode coupling could be 
due to electron-phonon or electron-electron interac- 
tions. Electron-phonon interaction, which could be en- 
hanced due to the proximity of the transition energy 
(250 meV) to that of the zone-center optical phonon 
(196 meV), must be excluded due to the different 
symmetry (this phonon is Raman active). Splitting 
due to Coulomb interaction can be evaluated to be 
0.04(e 2 /47T£o^-F)(7i73/7o)r^ i- e - omv about 3 meV in 
the absence of dielectric screening, e 2 / 'AtteqUvf = 2.2. 
Hence, the mechanism of the possible strong coupling 
between the L_i(_2) — > L-2(i) and L (_4) — > L 4 ( ) transi- 
tions is a puzzle which remains to be clarified. 



IV. CONCLUSIONS 

We probed graphene bilayers embedded in multilayer 
epitaxial graphene on the C-terminated surface of SiC. 
These inclusions can be viewed as AB-stacked faults in an 
otherwise rotationally-stacked multilayer graphene struc- 
ture and enable spectroscopic studies of unperturbed 
graphene bilayers. The "electronic quality" of these bi- 
layers is comparable or even better than that of the bi- 
layers obtained by exfoliation or by epitaxial growth on 
the Si-terminated surface of SiC . 34 ' 35 This way, we could 
trace the inter-band inter-LL transitions in the graphene 
bilayer for the first time, and thus supply data com- 
plementary to the cyclotron resonance absorption (i.e., 
intra-band inter-LL transitions) measured on the exfoli- 
ated bilayer by Henriksen et al^ We could also clearly 
visualize the departure of Landau levels in the graphene 
bilayer from the linearity in B, which clearly sets limits 
for the parabolic approximation of electronic bands in 
this material. 
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